Zero-mode analysis of quantum statistical physics 
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We present a unified formulation for quantum statistical physics based on the representation of 
the density matrix as a functional integral. We identify the stochastic variable of the effective 
statistical theory that we derive as a boundary configuration and a zero mode relevant to the 
discussion of infrared physics. We illustrate our formulation by computing the partition function 
of an interacting one-dimensional quantum mechanical system at finite temperature from the path- 
integral representation for the density matrix. The method of calculation provides an alternative 
to the usual sum over periodic trajectories: it sums over paths with coincident endpoints, and 
includes non- vanishing boundary terms. An appropriately modified expansion into Matsubara modes 
provides a natural separation of the zero-mode physics. This feature may be useful in the treatment 
of infrared divergences that plague the perturbative approach in thermal field theory. 



I. INTRODUCTION 

Quantum statistical physics provides the conceptual 
and computational framework for the treatment of inter- 
acting many-body systems in thermal equilibrium with 
a heat bath. It describes a great variety of systems and 
phenomena over a wide range of scales, thus covering 
practically all areas of physics, from cosmology to par- 
ticle physics, with an extensive number of examples in 
condensed matter. 

Its broad spectrum of applications includes special rel- 
ativistic systems. For those, the quantum statistical 
treatment is known as finite-tcmpcrature field theory [1] , 
which is synonymous to quantum statistical field the- 
ory. In fact, just as quantum statistical mechanics adds 
stochastic probabilities to the quantum probabilities of 
quantum mechanics, finite temperature field theory does 
likewise with respect to quantum field theory, the natural 
combination of special relativity and quantum mechan- 
ics. 

Quantum statistical physics can be formulated in the 
language of euclidean (imaginary-time) functional in- 
tegrals, quite natural for finite-temperature field theo- 
ries, but also applicable in nonrelativistic contexts, since 
quantum mechanics can be viewed as a field theory in 
zero spatial dimensions. Functional integrals thus furnish 
a powerful unifying formalism to compute correlations. 

Indeed, the formalism expresses the partition func- 
tion of a given system as a generating functional of eu- 
clidean Green functions (correlations), similarly to zero- 
temperature field theory, but with the time direction 
made compact, and specific boundary conditions that 
constrain the domain of field configurations in the func- 



* abessa@if.usp.br 
t aragao@if . ufr j . br 
^ fraga@if . ufr j . br 



tional integral. 

In order to best exploit the unifying aspect of the for- 
malism, we propose to view the partition function as an 
integral of the diagonal density matrix element of the the- 
ory. The integral is performed over the stochastic vari- 
able which characterizes the representation of the matrix 
element (position representation for quantum mechanics, 
field representation for second quantized quantum me- 
chanics or field theory). The density matrix is just the 
Boltzmann operator, whose Hamiltonian operator may 
cither come from (second quantized) quantum mechan- 
ics or from field theory 

The density matrix element may be written as a func- 
tional integral in euclidean time t which starts from a 
given configuration (a point in quantum mechanics, a 
field configuration in field theory) at r = 0, and returns 
to that same configuration at r = (3h. Clearly, there is a 
sum to be performed over configurations which coincide 
at the endpoints of the euclidean time interval. Those 
boundary configurations are identified as the stochastic 
variable of the remaining integral. Furthermore, it will 
be shown that they are zero modes of a modified Mat- 
subara series. The theory of the zero modes is the effec- 
tive stochastic theory obtained from the original Hamil- 
tonian. 

This density matrix method for deriving an effective 
stochastic problem was already used in [2], where it led 
to the construction of dimensionally reduced effective ac- 
tions. Likewise, in Ref. [3] we used the density matrix 
method, and a scmiclassical approximation, to investi- 
gate the thermodynamics of scalar fields. In the latter 
reference, the boundary configuration of the field (writ- 
ten as a constant plus gaussian fluctuations) played a 
very nontrivial role in the semiclassical formulation for 
the partition function of scalar fields. 

Our prior uses of the method drew upon generaliza- 
tions of a thermal scmiclassical treatment previously de- 
veloped and applied to quantum statistical mechanics, 
which produced excellent results for the ground state en- 
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ergy and the specific heat of the anharmonic (quartic) 
oscillator [4-6]. In those prior uses, however, the connec- 
tion with the stochastic variable was indirect, through 
classical solutions of equations of motion. Now, we have 
opted instead for a direct connection with the stochas- 
tic variable, to be identified with a zero mode, because 
we believe this is physically relevant to the discussion 
of infrared problems in several applications of quantum 
statistical physics. 

In the context of finite-temperature field theory, for 
instance, it is known that a naive implementation of per- 
turbation theory for the calculation of Feynman diagrams 
is ill-defined in the presence of massless bosons. This 
is due to the appearance of severe infrared divergences, 
brought about by the vanishing bosonic Matsubara mode 
in thermal propagators. The divergences plague the en- 
tire series, making it essentially meaningless [1, 7, 8]. 

As a result, one is forced to resort to resummation 
techniques that reorganize the perturbative series, and 
resum certain classes of diagrams, in order to extract 
sensible results. There are several ways of performing 
resummations, and rewriting the degrees of freedom more 
efficiently in terms of quasiparticles; we refer the reader 
to the reviews [7, 8] , and to Ref . [3] , for a discussion and 
a list of specific references. 

All such techniques are designed to partially tame the 
infrared divergences, creating a non-zero domain of va- 
lidity for weak-coupling expansions^. Nevertheless, the 
zero-mode problem remains, and the region of validity 
of resummed perturbation treatments can not be indefi- 
nitely enlarged. 

The infrared problems of finite temperature field the- 
ory resemble those encountered in the functional inte- 
gral treatment of various condensed matter systems [10], 
where they are connected to collective excitations. This 
suggests that one should explore the unifying feature of 
the formalism, which is its statistical nature, to profit 
from parallel physical interpretations and insights. This 
is yet another reason for resorting to the density matrix 
formulation. 

In the present article, we resume the study of quan- 
tum statistical mechanics, viewed as an exercise in zero- 
dimensional finite-temperature field theory. Our objec- 
tive is to show that boundary configurations are the zero- 
mode stochastic variables of the effective statistical the- 
ory, and to compute that theory in various approximate 
schemes. Although we restrict our analysis to quantum 
mechanical examples, this should be regarded as a pre- 
liminary to the field theory case. 

We compute the partition function of an interacting 
one-dimensional quantum mechanical system at finite 
temperature from a path-integral representation for the 
density matrix. The method of calculation that we pro- 



pose provides an alternative to the usual sum over pe- 
riodic trajectories: it sums over paths with coincident 
endpoints, and includes the contribution of non- vanishing 
boundary terms. Indeed, an appropriately modified ex- 
pansion into Matsubara modes provides a natural separa- 
tion of the zero-mode physics, which is connected to the 
infrared problems in thermal field theory and elsewhere, 
and relates the zero mode to the boundary value of the 
quantum- mechanical coordinate. 

The paper is organized as follows: in Section II, we 
present the proposed modified series expansion in Mat- 
subara frequencies, and make a detailed comparison of 
our density matrix procedure with the standard one, for 
clarity; in Section III, we relate the boundary value of the 
quantum-mechanical coordinate (the stochastic variable) 
to the zero mode; In Section IV, we compute several ther- 
modynamic quantities for the quadratic and the quartic 
potentials, and compare our results to the semiclassical 
findings of [4], and to some exact results; Section V con- 
tains our conclusions and outlook. 



II. A MODIFIED SERIES EXPANSION 

The partition function Z for a one-dimensional 
quantum-mechanical system in contact with a thermal 
reservoir at temperature T = \/{kBl3) may be written as 



(1) 



where p(/3;xo,xo) is the diagonal element of the density 
matrix. In a path-integral formulation, the density ma- 
trix p is obtained from the imaginary-time evolution of 
the trajectory x : [0, (31t\ M determined by an euclidean 
action S: 



p{f3; xo,xq) 



[Vx] 



-s[x]/h 



x{0}=x{l3h)=xo 



S[x\ = dT 
Jo 



i ■mx^{T) + V{x{t)) 



(2) 



(3) 



where V is the potential^. In such a formalism, it is 
natural to regard xq as the effective stochastic degree of 
freedom of the theory, and to compute the distribution 
of the variable xo by integrating over the auxiliary vari- 
able x, with boundary conditions determined by Xq. In 
this article, we want to show that the effective theory for 
the static variable xq not only underscores its statisti- 
cal mechanical character, but can also be identified with 
the theory of the zero mode of an appropriately modified 
Matsubara expansion. 



^ In the case of thermal QCD, for instance, the domain of validity 
of the naive perturbative expansion is the empty set [9] . 



^ In order to simplify the notation, we will adopt natural units 
where h = 1 and ks = 1 . 
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We can show that the formalism requires a modified 
Matsubara expansion for x by calculating the partition 
function of the harmonic oscillator in two different man- 
ners. The action S, in this simple case, is given by 



S[x] 



(4) 



Without any loss of generality, one can take m = 1 in the 
previous equation by rescaling x by a factor 1 / \/rn. One 
may compute the partition function using the following 
path integral: 



J [Dx] e 



(5) 



The condition a;(0) = a;(/3) is usually implemented by 
expanding the path x in Matsubara modes: 



x{t) 



E 



Cnn 6 



with CiJri 



n— — CX3 



27rn 



(6) 



The standard procedure is, then, to integrate Eq. (4) by 
parts: 



S{x\ 





r d^ ,1 


/ dTx{T) 

Jo 





xir) 



x{t)x(t) 



(7) 



Using the periodicity of Eq. (6), one shows that the 
boundary term in (7) is zero. The remaining term can 
be cast in the form 



S\x\ 



drdr' x{t)K{t, T')a;(T'), 



(8) 



where the Green function K{t,t') ^ Ap{t ~ t') is such 
that 

-^+c.2 Af{t) = S{t), (9a) 
Af{t ^ (3) ^ Af{t) . (9b) 
This corresponds to the standard free propagator [1] 



(10) 



where n(Lo) is the Bose-Einstein distribution. Using Ap, 
it is straightforward to obtain the partition function. Up 
to an infinite constant, one obtains the well known result 
for the free energy 



1 u 1 , 



(11) 



There exists, however, an alternative path-integral 
method to calculate the partition function of the har- 
monic oscillator, which has the advantage of being free 



of spurious divergences. We start with (1), and decom- 
pose the configuration x as 



x(r) = Xc{t) +y{T) , 



(12) 



where Xc is the solution of the following Euler-Lagrange 
equation: 



Xc{0) = Xcif3) ^ Xo 



which is given by 



Xcir) = Xo 



cosh[u;(T-/3/2)] 
cosh([j/3/2) 



(13a) 
(13b) 

(14) 



It is easy to show that 



1 



S[x] ^ -xo[xi(3) ~ m] + S[y] 



XQUj{cosh f3uj — 1) 
sinh Plu 



S[y] 



(15) 



Using that y vanishes at r = 0, /3, one rewrites the action 
of the trajectory y as 



1 



S[y]^:: / drdr'y{T)G-\r,T')y{r') 



(16) 



where the Green function Gy satisfies the following equa- 
tion: 



dr2 



Gy{r,r') = S{t^t'), 

Gy{0,T')^Gy{l3,T')=0. 

One can show [4] that 

sinh[ti>(T> — /?)] sinh(/3r<) 



G,(r,r') = 



Lu sinh(w/3) 



and 



'^''^^y^ 2n sinh(c^/3) 



(17a) 
(17b) 

(18) 

(19) 



In Eq. (18), T> = max(r, r') and t< = min(T, t'). The 
path integral over y produces, then, 

Z = {detGy^^^ J da;ge-"[('=°"'''''^-i)/'''"'''^"l^'' (20) 

^ {2 cosh Plu -2)-^/^ . (21) 

Using (2cosha;/3 — 2) = [1 — exp(— /3a;)]^ exp(/3w). and 
taking the logarithm we recover Eq. (11). 

Notice, however, that the two methods treat in a dif- 
ferent way the quantity 



x{P) - i(0) 



(22) 
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where 



FIG. 1: Schematic construction of the extension of a config- 
uration X defined originally in [0, /3] to the interval [— /3] . 
Notice the odd character of the extension of a; — a;o. 



Let us focus on the usual Matsubara expansion of the 
classical path (Eq. (14)). Being periodic, the correspond- 
ing series has a periodic derivative, so that the quan- 
tity (22) is put to zero. That was used to eliminate the 
boundary term in Eq. (7). However, we see from (15) 
that x'c(/3) — a;c(0) ^ 0, otherwise one would obtain an 
incorrect result for Z. 

The second method of calculation sums over paths 
which have coincident endpoints, and thus includes paths 
that are not periodic. It is easy to realize that we in- 
troduce an error in the value of the action of any non- 
periodic path X by evaluating the action of the associate 
Matsubara series as given by Eq. (6) . The reason is that 
the usual Matsubara expansion is not suitable to handle 
boundary terms involving the derivative of a; at and /3. 
We stress that, as long as the only condition over paths 
in the path integral is to have coincident values at and 

non-periodicity is allowed. 

One can remedy this situation by defining an extension 
of X to [— /3,/3] as follows (see Fig. 1): 



x{t) = x{t) X[o,/3] + (2x0 - x{-t)) X[-/3,o) 



(23) 



where xa denotes the characteristic function of A C M. 
Notice that i(— r) — xq = —{x{t) — xq). The odd char- 
acter of the extension implies that x — xq has a series of 
sines: 



Xn smw„ T . 



with LjJ„ 



rnr 



(24) 



The notation with a dot is a reminder that the r.h.s. of 
Eq. (24) is a Fourier approximation of the l.h.s. The 
desired restriction to [0, /3] is: 



x{t) = xq + Xn sina)„ 



(25) 
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dr x{t) sina)„T . 



(26) 



The series (25) does not impose conditions on the deriva- 
tive of X. It is easy to convince oneself that evaluating 
the euclidean action S with the r.h.s. of Eq. (25) one ob- 
tains exactly S[x]. In particular, we have for the classical 
solution: 



Xcir) = Xq - Xq 



n odd 



SmWnT , 



mr {uj^^ + up') 
so that the quantity XcXfi) — a;c(0) is given by 



(27) 



E—AxQu'^ujn , . „ ^ 2a;ocj(cosh/3a; — 1) 
^ . 2 I . .2^ (cosw„/3 - 1) = 



sinh Puj 



(28) 



as implied by Eq. (15). Obviously, the series defined 
by Eq. (25) does not correspond to the most general 
trajectory in Z. In particular, that series imposes a seri- 
ous restriction on the second derivative of x, by forcing 
it to vanish at and /3. However, at least for the cal- 
culation of usual action functionals, the boundary terms 
do not involve second derivatives and the aforementioned 
restriction is actually unimportant. 



III. THE BOUNDARY VALUE AS THE ZERO 
MODE 

The first remarkable property of the decomposition 
(25) is that the boundary value xq is the zero, static com- 
ponent of X which, in the context of thermal field theo- 
ries, plays a special role in the infrared (low-momentum) 
limit of the theory. As mentioned previously, the in- 
frared physics of bosonic theories is not accessible to plain 
perturbation theory, and resummation techniques are re- 
quired in order to produce sensible results. In common, 
all such techniques render a special treatment to the zero 
mode'^. 

It is important to stress that the zero mode xq is not 
the zero mode of the usual Matsubara expansion. In 
fact, from Eq. (6), one can notice that xq involves the 
coefficients of all Matsubara modes: 



Xq 



E 



n— — oo 



(29) 



The present identification of the zero mode with the 
boundary field in the context of quantum statistical me- 
chanics aggregates to the zero-mode physics yet another 



^ Rigorously speaking, xq is not a true Fourier mode since it is not 
orthogonal to sintD„T in [0,/3], but in [—/3,l3]. With this caveat 
in mind, wo will continue to refer to xq as a mode of x. 
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important feature: the connection with the classical be- 
havior at high temperatures. That regime is character- 
ized by the condition 



(30) 



(or /3a; ^ 1, in natural units), where hu) is the typical 
spacing between quantum levels. In order to see why the 
static mode governs the classical limit of the partition 
function, let us consider Xc, the solution of the equation 
of motion: 



(fxc 



xc{0) = Xci/3) 



y'(x,(r)) -0. 

xo . 



(31a) 
(31b) 



In the high-temperature regime (/3a; ^1), thermal fluc- 
tuations dominate. In this limit, Xc — > xq and S{xc) — > 
(3V{xo), and we obtain for the partition function 



lim Z 

T->oo 



/ dxo e 



(32) 



where is a normalization factor which incorporates 
quantum fluctuations. A proper normalization of Z leads 
to iV^ = mkBT/{2TTh'^), rendering the semiclassical ex- 
pression 



lim Z = 



dp 
2Trh 



dxo e 



-H/kBT 



(33) 



with H = /2m + V{x{)), m being the particle mass. 

As discussed in [11], for /3a) ^ 1 the correlation {x^ij)) 
follows the classical linear scaling with T. This behav- 
ior is entirely associated with the static component of 
x, and represents a problem for the high-temperature 
limit of the usual perturbative expansion. In contrast, 
the subtracted thermal propagator (without the contri- 
bution from the static mode) goes to zero as T increases, 
so that one should strongly improve the convergence of 
the perturbation expansion by calculating diagrams with 
that subtracted propagator. 

We conclude that, either to control infrared diver- 
gences in thermal field theories or to describe the classical 
limit in quantum statistical mechanics, we need a sepa- 
rate treatment of the zero mode. The main advantage of 
our approach is to provide this separation in a natural 
fashion. 

In order to explore the zero-mode physics, it is conve- 
nient to write 



x{t) = .To + yir) , 



(34) 



where y corresponds to the sinusoidal modes of Eq. (25) 
and has the following property: 

2/(0) = 2/(/3) = . (35) 

In terms of y one can express Z conveniently as 



dxQ 



[Vy] e 



-S[xo+v] 



(36) 



and different strategies can be used to handle S[xq + y] 
in order to obtain an effective theory** for a;o: 

dxQ 



Z = / • " e-^ . (37) 

J ^2T:h? /{mkBT) 

Notice that the intermediate step, including the paths y, 
is built using quantities which vanish at the boundary 
of imaginary time, and the contribution from problem- 
atic static quantities is factorized in the final integration 
over xq. Besides, for interacting systems, the interaction 
potential evaluated at xir) = xq + y^r) can produce a 
2/-dependent quadratic term for the zero mode xq, which 
could prevent possible divergences in the effective theory. 

In the sequel, we use a very simple expansion around 
xq in two cases: the harmonic potential, as a consistency 
check, and the single-well quartic potential, in order to 
test the quality of our alternative approach and compare 
it to previous descriptions. 



IV. APPLICATIONS AND COMPARISON TO 
PREVIOUS RESULTS 

A. The quadratic case revisited 

The action is again the one given by Eq. (4). Now, we 
use the decomposition of x given by Eq. (34) to obtain 



S[x] 



drdr' {yG-^y + 2uj^xoy} 



2 + ^ f dT{{y + u;^xoGy) Gy^x 



{y + uj^xo Gy) } 



4^2 



UJ X: 



drdr'Gy , (38) 



where Gy was defined in (17). A simple calculation gives 

^ ^ 'r^ / M ^ 2(cosh/3^-l) 
drdr Gy[T, r ) = — ^ . . „ . (39) 

Shifting the variable of integration to 

C(r) = y{T) + cj'xo I dr' Gy{T, r') , (40) 



which obeys QiQ) = CiP) ~ 0- a-nd performing the gaus- 
sian integration, we obtain 

Z = {det Gy f^^ [ rfa;ge-'^('=°'^'^^'^-i)/^'"'>'3'^)^o . (41) 



Integrating over xq and taking the logarithm, one repro- 
duces Eq. (20). Notice that the propagator used in the 
perturbation expansion of the integral over the path y is 
well-behaved in the limit of small u: 

|G,(r,r')| < ^^^^^ ^. (42) 



y(o)=y(/5)=o 



* With the normahzation factor used in (37), one has the property; 
Voff(xo) V{xo) when T — > 00, as follows from (33). 
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B. The anharmonic oscillator 



Let us now consider the following action: 



S[a 



dr 



-mx (r) H — muj x (r) H — a; (r) 
2 2 ^ ' 4 ^ ' 



(43) 

It is convenient to define the dimensionless variables: q = 
y^\/{mu;'^)x, 6 = uiT, e = ljP and g = \/{m^uj^). The 
associated action for the variable q is: 



1 



I 



I 



I 



S[q] = -g [l^i'iO) + 2«'(^) + j • (44) 

Using the decomposition q{0) = go + v^^(^)' write 
for the action: 



S[q] = S2[qo,v] +Si[qn,ii\ , 



(45) 



where ^2 contains 5[(jo], the kinetic term, and those 
which are linear or quadratic in rj. The remaining contri- 
butions to S[q] are collected in 5/. It is a simple matter 
to show that 



de^ ^ 

2a 

^71 



m \ , (46) 



and 



Si[qa,'ri\ = 



d6 



(47) 



where = Sgg + 1 and a = qo + qo- One obtains 
a quadratic approximation for the partition function by 
considering 



dqo J [Pr/le-^^I*^"! 
>)(o)=»)(e)=o 



(48) 



Proceeding as in the previous section, we perform the 
gaussian integration over the variable 



(49) 



m^rm + ^ / d9'G,{9,9') 

Va Jo 

where the Green function G,,, defined by 



de^ 



G,,i9,9') = S{9-9'), 



G^{0,9')^G^{e,9')^0 



is given by 



Gr,{9,9') 



sinh[tZJ(6'> - e)] sinh(LZJ6'<) 
tUsinh(ZI) 0) 



(50a) 
(50b) 

(51) 



The result for is then 



Z2 = J dqo (dot G^)i/' exp{-5[go] + S J , (52) 



where 



and 



det G„ = r-r-, — ^ 

' 27r sinh(ZZ;e) 



S„ = — / d9d9'GJ9,9') 



(53) 



(54) 



To go beyond the quadratic order in 77, we use the fol- 
lowing approximation: 



Jo 



V5 W(^) + ^^'(^?) • (55) 



Therefore, Z w Z2 + S^^^Z, where 

6^'^z = - f dqo r d9 [jdqoin'm + \{v\e)) 



(56) 



with 



(77(^1 )... 77(0fc )) = j [I?r;]e-^=l*'"l77(ei)...7?(0fe) . 

»,(0)=r,(e)=0 

(57) 

As usual, we can introduce a linear coupling of 77 with an 
external current j: 

Z2[j]^ J dqo J [p77]cxp|-52[go,?7] + J d9j{9U9) 



r,(0)=r,(e)=0 



(58) 

and obtain the correlations (57) as functional derivatives 
of Z2[j] with respect to j. With minor modifications on 
the recipe for Z2, one obtains 

Z2[:j] = J dqo (det G,)'/' cxp{-5[go] + Srjj]} , (59) 



where 



1 



^„[j] = -J^ d9d9'a{9)Gr,{9,9')a{9') , (60) 
and (j{9) = a/^/g — j(9). Finally, we use 



6f{9) 



cS„b1 ^ /3(0)+3/o(0)G' 



3=0 



(61) 



and 



e^.W =/4(0)-)-6 70^(0) G,(0,0) 

+ 3Gl{9,9), (62) 
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FIG. 2: Plot of r,„i„ as a function of the coupling constant A 
(m = l,ui = 1). The quadratic approximation is applicable 
for temperatures larger than T„i„. 

where /o(^') plays the role of an expectation value of rj{9): 

io{0)^-^ I de'Gr,{e,d') 

Vff Jo 

a cosh-UjO — smhioO tanhLj0/2 — 1 
= =2 ■ (63) 

The remaining integrations over 9 can be done analyti- 
cally, providing all the ingredients for the calculation of 
the first correction to Z2. 

Notice that the integrand of (52) is dominated by the 
vicinity of (?o = 0, where the quadratic action reaches its 
minimum value. For small temperatures, that integrand 
becomes sharply peaked at zero, while for large T a broad 
range of go substantially contributes. Therefore, one can 
estimate the importance of the correction term (56) using 
the following quantity: 

(5z(0,,))=|£rf^G^(^,^) -^^J^. (64) 

The temperature T,,,;,, where (5/(0,?/)) goes to 1 is dis- 
played in Fig. 2 as a function of the coupling constant 
A. r„i„ provides a rough estimate for the validity of the 
quadratic approximation. From Fig. 2, we see that T^^^ 
is of order 1 for A as large as 50. 

The results for the free energy and for the specific heat 
are displayed in Fig. 3 and Fig. 4, respectively. In the fig- 
ures, classical corresponds to calculation using Eq. (33), 
gitadrafic uses Eq. (52). improved includes the correction 
corresponding to Eq. (56), semiclassical refers to the 
scmiclassical result obtained in Ref. [4], 1-loop stands 
for the usual 1-loop perturbative result, and exact cor- 
responds to results from Refs. [12, 13] combined with 
WKB estimates. 

The first remarkable thing in the plots is the failure 
of the perturbative result in the high-temperature region 
(T 3> m). As seen in Ref. [4], that region is well de- 
scribed by the semiclassical curve, and asymptotically, 
by the classical one. This is evident from both plots. It 
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FIG. 3: Free energy of the anharmonic oscillator as a function 
of the temperature for A = 0.4, m = 1 and to = 1 (see text 
for details). 



is interesting that, for g = 0.4, the agreement extends 
to temperatures down to T « m/2, a region where the 
character of the system is far from classical. However, it 
is surprising that our simple calculation using (52) works 
as well as the semiclassical one in that region. 

Indeed, the cited semiclassical calculation is based on 
a quadratic expansion around exact solutions of (31) for 
the quartic potential. In practice, one has to deal with 
rather involved Jacobi elliptic functions [4, 14]. In con- 
trast, using the expansion (34) no knowledge of classi- 
cal solutions is demanded. This simplification represents 
an economic alternative which can be crucial in contexts 
where exact classical solutions are not available. From 
the plots, one also concludes that, in the region where the 
approximation is supposed to be good, the improved cal- 
culation exhibits a stronger convergence, indicating that 
the approximation is consistent. 

A detailed look at the quadratic curve in the free en- 
ergy plot (Fig. 3) reveals a non-monotonic behavior for 
low temperatures: the free energy passes by a maximum 
and decreases towards the free value 0.5. This is a gen- 
eral feature regardless of the value of m and A, and it is a 
clear signal of the breakdown of the approximation below 
T,-„i„ (see Fig. 2). The improved free energy diverges at 
T,-„i„, because at this value the sum Z2 + vanishes. 
The exact curve for the free energy reaches its maximum 
value at the ground state energy, where it rests down 
to r = 0. Using the maximum value assumed by the 
quadratic curve as an approximate value for the ground 
state energy, we obtain unexpectedly reasonable results, 
as shown in Table I. 

Fig. 4 compares different predictions for the specific 
heat for A = 0.4, to = 1 and u) = 1. Notice that all 
curves, except for the perturbative one, have the cor- 
rect high-temperature limit. The anomalous behavior 
of the quadratic approximation near T„i„ is even more 
evident in this plot. The improvement obtained in the 
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FIG. 5: Specific heat of the anharmonic oscillator with lj = 0, 
m = 1 and different values of A. 



low-temperature regime when one corrects the quadratic 
calculation with (56) is also clearly shown. 

Finally, Fig. 5 displays the quadratic result for the 
specific heat when A = 0.4 and m = 1 in the limit of zero 
u}. For large temperatures, the curve converges to the 
classical value 0.75. As discussed in Section II, with the 
separation of the zero mode, the finiteness of the theory 
for the 77 variable is not affected in the infrared limit, and 
we see in this specific case that the interaction makes the 
effective zero-mode theory well defined. 



V. CONCLUSIONS AND OUTLOOK 

The theoretical description of the thermodynamics of 
a given interacting many-body system kept in thermal 
equilibrium by the contact with a heat bath is an issue 
of capital importance in most realms of physics, from 
cosmology to condensed matter settings. A very conve- 
nient framework is provided by finite-temperature field 
theory, where one can use the imaginary-time formalism 



A 


£0 (exact) 


iJo (quadratic) 


Error(%) 


0.008 


0.501 


0.505 


0.8 


0.04 


0.507 


0.518 


2.2 


0.4 


0.559 


0.584 


4.5 


1.2 


0.638 


0.662 


3.8 


2.0 


0.696 


0.718 


3.2 


4.0 


0.804 


0.818 


1.7 


8.0 


0.952 


0.958 


0.6 


200.0 


2.500 


2.450 


2.0 



TABLE I: Ground state energy of the quartic oscillator for 
different values of the coupling A (m = a; = 1). Exact data 
were extracted from Ref. [13], whereas iSo (quadratic) was 
obtained using Eq. (52). 



and make a direct connection to the powerful formulation 
of thermal averages in terms of functional integrals. Nev- 
ertheless, whenever collective (masslcss) bosonic modes 
arc relevant, as is the case in practically any plasma, 
infrared divergences are unavoidable in a plain pertur- 
bative treatment, forcing a reorganization of the series 
that grants special attention to the zeroth Matsubara 
mode. As discussed in this paper, we also need a sepa- 
rate treatment of the static mode already in the context 
of quantum statistical mechanics in order to improve the 
convergence of the perturbativc expansion in the high- 
temperature limit. 

The choice of the traditional Matsubara mode expan- 
sion, combined with the trace requirement that the paths 
(configurations) coincide at the two boundaries of the 
compact euclidean time, implies a restriction to periodic 
paths in the process of functional integration to compute 
the partition function. However, rigorously speaking, the 
paths need not be periodic in the path-integral represen- 
tation of the partition function. Indeed, the necessity of 
paths which are not periodic is particularly clear in the 
formulation of the partition function as an integral of the 
diagonal density matrix clement of the theory. We have 
shown that, in order to allow for non-periodicity, one has 
to define a modified expansion into Matsubara modes. 

In this paper, we have followed our alternative strategy 
in the case of quantum statistical mechanics, viewed not 
only as a toy model for the case of finite-temperature 
field theory, but also as a prototype for various rele- 
vant systems in statistical physics. More specifically, we 
have explored an alternative way of computing the parti- 
tion function in the path-integral formalism that includes 
non-periodic trajectories, for which the usual Matsubara 
expansion leads to an incorrect result for the action of the 
trajectory. As was shown above, one can properly incor- 
porate the contribution of non-periodic paths by using a 
modified Matsubara series expansion in which the zero 
mode is identified with the boundary value of the path. 
The latter turns out to be the stochastic variable that 
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survives in the final effective theory. 

The approach proposed in this paper has, thus, the ad- 
vantage of providing a natural separation of the physics 
of the zero mode. In fact, we built a very simple effec- 
tive theory for the zero mode in the nontrivial problem of 
the anharmonic oscillator, obtaining very precise results 
for practically the whole range of temperatures, and in 
a framework that is free of spurious divergences. In par- 
ticular, we were able to describe the semiclassical (high- 
temperature) limit of the anharmonic oscillator without 
any information about the classical solutions. We expect 
that other potentials can profit from the simplicity of the 
present method in quantum mechanical systems. 

In the case of finite-temperature field theory there are, 
of course, subtle issues of renormalization, which are ab- 



sent in the quantum mechanical setting, that will have 
to be addressed. Nevertheless, the results obtained in 
quantum statistical mechanics are encouraging, and we 
hope that this new perspective may shed light onto the 
problem of infrared divergences. Results in this direction 
will be presented elsewhere [15] 
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